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Abstract. We calculate the optical depth and the number of 
events due to gravitational microlensing towards the Galactic 
bulge, the spiral arm directions 7 Scutum, (3 Scutum, 7 Nor- 
mae, d Muscae and some dwarf galaxies in the halo of the 
Galaxy. 

Using the events found by the MACHO collaboration dur- 
ing their first year of observation towards Baade's Window we 
estimate the mass functions for the bulge and disk populations 
following the mass moment method. We find that the mass 
function can be described by a decreasing power-law with slope 
a ~ 2.0 in both cases and a minimal mass of ~ 0.01 Mq for 
the bulge and ~ 0.02 M Q for the disk, respectively. Assum- 
ing that the obtained mass function for the disk is also valid 
in the spiral arms, we find that the expected number of events 
towards the spiral arms is in reasonable agreement with the ob- 
servations. However, the small number of observed events does 
not yet constrain much the different parameters entering in the 
computation of the mass function. 

To study the influence of the Magellanic Clouds on the 
shape and the velocity dispersion in the halo we perform a N- 
body simulation. We find that their presence induces a slight 
flattening of the halo (qn — 0.8). As a result the expected num- 
ber of microlensing events towards some targets in the halo, 
such as the LMC or the SMC, decreases by about 20%, whereas 
due to the the modification induced on the velocity dispersion 
the event duration increases. 

Key words: dark matter - Galaxy: stellar content, structure - 
microlensing - stars: low-mass, brown dwarfs 



1. Introduction 

During the last ten years amazing progress was made by ex- 
ploring the dark component of our own Galaxy by means of 
gravitational microlensing, as proposed by Paczyfiski in 1986. 
Microlensing allows the detection of MACHOs (Massive As- 
trophysical Compact Halo Objects) in the mass range 10~ 7 < 
M/Mq < 1 (De Rujula et al. |1992|) located in the Galactic 



1991,Griestetal. 1991). Microlensing searches have nowadays 



become an important branch in astrophysics, especially for the 
study of the structure of the Milky Way and even for globular 
clusters (Jetzer et al. 1998). 

Today, more than half a dozen groups are active in observ- 
ing microlensing events towards different lines of sight, and 
they reported so far several hundreds of events, most of them 
towards the Galactic bulge (Alcock et al. 1997a, Udalski et al. 



1994 , Ala rd et al . 1997 ), so me tow ards the spiral arms (Al- 
cock et al. |l997aj Deru e et al . |l999j), about 15 events towards 



halo, as well as in the disk or bulge of our Galaxy (Pazcyhski 
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the LMC (Alcock et al. [199J , |1997bj Aubur g et al |1993| ) and 
two events toward s the SMC (Alcock et al. 1997c , Palanque- 
Delabrouille et al. 1997 , Udalski et al. 1997 ). The necessity to 
extend the microlensing target regions, as for instance towards 
M31 (Crotts & Tomaney 1996, Ansari et al. 1999), is widely 
recognized in order to better characterize the MACHO distri- 
bution throughout the Milky Way and the Galactic halo. 

In this paper we calculate the optical depth and the number 
of events due to gravitational microlensing towards the bulge, 
the spiral arm directions 7 Set, j3 Set, 7 Nor, d Mus and some 
dwarf galaxies in the halo of the Galaxy. To that purpose we use 
the me thod o f mass and time moments (De Rujula, Jetzer and 
Masso 1991 ). Based on a model of our Galaxy and taking into 
account the different populations of objects that can act either 
as lenses or as sources in microlensing events, we compute the 
different mass and time moments, which are needed for the 
determination of the various microlensing quantities. Due to 
the extension of the bulge we include in the calculation of the 
optical depth and the expected number of events the fact that 
the number density of sources varies with distance. 

The 41 events published by the MACHO collaboration to- 
wards Baade's Window allow us to calculate the mass functions 
for the objects acting as lenses towards the Galactic bulge. The 
result is described by a decreasing power-law with M m i n ~ 
0.012M© for lenses in the bulge and M mm ~ 0.021M Q for 
lenses in the disk and a slope a ~ 2.0 in both cases. Moreover, 
assuming that the disk mass function is the same also in the spi- 
ral arms, we find that the calculated number of events towards 
the spiral arms is in good agreement with the events reported 
by the EROS II collaboration. 



2 



Microlensing Events and Dark Objects by Mass and Time Moments 



However, the small number of observed events towards the 
spiral arms does not yet constrain much the different parame- 
ters of the mass function. 

We emphasize that, due to the small number of events at 
disposal, the mass functions we found for the bulge and the 
disk should be considered as an illustrative example of how the 
mass moment method can be used to get useful information 
on the physical parameters. Indeed, when the many new events 
which have been observed in the meantime will be published, 
more accurate results will be obtained. 

To study the influence of the Magellanic Clouds on the 
shape and the velocity dispersion in the halo we perform a N- 
body simulation. We find that their presence induces a slight 
flattening of the halo (qn — 0.8) starting from an initially 
isotropic distribution. As a result the expected number and the 
duration of microlensing events towards targets in the halo gets 
modified. 

The paper is organized as follows: in Sect. 2 we present 
our Galaxy model, in Sect. 3 we introduce the method of mass 
and time moments. The calculation of the various microlens- 
ing quantities towards the Galactic center and the spiral arms 
is done in Sect. 4, while Sect. 5 is dedicated to different targets 
in the halo of the Galaxy. A short discussion in Sect. 6 con- 
cludes the paper. The N-body simulation of the Galactic halo is 
presented in the Appendix. 
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where R, z are cylindrical coordinates. Here we adopt i?thin — 
0.3 kpc for the thin and -f/thick — 1 kpc for the thick disk com- 
ponent, whereas h ~ 3.5 kpc for the length scale of the two 



disks. Following Gates et al. (1995 1, we take for the local sur- 
face densities Ethin — 25 M Q /pc 2 and Schick — 35 M /pc 2 . 
Rq = 8.5 kpc is the distance of the solar system from the Galac- 
tic center. 

The halo is assumed to be a self-gravitating isothermal 
sphere of an ideal gas in hydrostatic equilibrium. A slightly 
flattened halo can then be described by the density distribution 

PoH R c + H-a 
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where q B is the axis oblateness ratio (Binney and Tremaine 
R c ~ 5.6 kpc is the core radius and p n — 7.9 x 
~ 3 is the local dark mass density. 



1987| ), 

10~ 3 Mr. 



1q pc 

For the mass within 50 kpc we obtain for the bulge Mb = 
1.8 x 1O 1O M , disk M D = 5.2 x 10 10 M Q and halo M H = 
4.3 x 10 11 Mo, respectively. Thus the total mass is Mtot — 
5 x lO n M0 in agreement with e.g. Kochanek ( 1995 ). The cal- 
culated rotation curve agrees well with the measured values. 



2. Galaxy models 

Our Galaxy model consists of three components: a bulge, a disk 
and a halo, which we discuss in the following. 



2.1. Mass density 



Following Han and Gould (1995) we assume a triaxial bulge 
with a density law (Dwek et al. 1995) 1 



p B (x',y',z') = p oB e 8 /2 



M F 



6.57nabc 



-s 2 /2 



(1) 



with s 4 = (x' 2 1 a 2 + y' 2 /b 2 ^j +z' 4 /c 4 . M B ~ 1.8x 10 10 M Q 
is the estimated bulge mass. The length scales are: a = 1.58 
kpc, b — 0.62 kpc and c = 0.43 kpc (x' and y' are defined 
along and perpendicular to the bar shaped bulge in the Galactic 
plane). We take for the inclination angle between the bar major 
axis and the line of sight towards the Galactic center a value of 
$ = 20°. The influence of this parameter on the results will 
be discussed. 

We consider a disk model with a thin and a thick compo nent 
and a density d istribu tion as given by Bahcall et al. ( 1983 ) and 
Gilmore et al. ( |l989| ) 



p D (R,z) = - exp 



(R - Ro) 



x 



1 The integral exp( — s 2 /2)dx'dy'dz' evaluates to 6. 577ra6c. 
In the literature sometimes the wrong normalization l/(8nabc) is 
used instead. 



2.2. Mass distribution 

The mass density does not determine the MACHO number 
density as a function of mass alone. Assuming the mass- 
distribution to be independent of the position in the galaxy, the 
number density can be written as 



dn\ ^ / dn 
d PJi \ d P 



P°i 



dp,, 



(4) 



where /i is the MACHO mass in solar mass units. The subscript 
i stands for the bulge (B), the disk (D) or the halo (H). The MA- 
CHO number density per unit of mass, dn /dpi, is normalized 
as follows: 
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dn a \ 

— ) [idp = p 0i 
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The total mass distribution is just the sum over the components 
in Eq.(Q). We assume a maximal mass fj, max ~ 10 M Q for the 
stars. However, only the faint stars up to about fi up ~ 1 M 
can contribute to microlensing events. For the bright stars \i > 
1 Mq we assume a Salpeter IMF (dno/dp oc /i -2 35 ), whereas 
for the lenses we will either assume all objects to have the same 
mass p Qi and, therefore, the mass distribution is described by a 
delta function 



dn \ p Qi 
dp J i M Q fj, 0i 



S(p-p 0i ) 



or we assume a power-law 



(6) 
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Accordingly the factor C\ (oi) in the power-law is fixed by the 
normalisation condition (for i = B, D) 



(8) 



Assuming continuity for the mass function at jj, up — 1 M we 
get 



cl( ai ) 



Ma 



1.5809 



if a, £ 2 



(9) 



ln Mmm + 1.5809] 1 if a, = 2 



2.3. Velocity distribution 

For sources and lenses located in the bulge we assume that 
the velocity distribution along the various axes is Gaussian 
f{v y ,v z ) = f(v y )f(v z ), where 
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and a corresponding distribution for f(v z ). The mean velocities 
u a and v z are supposed to be zero and the dispersion velocities 
a y and a z are deduced from the tensor virial theorem (Binney 



& Tremaine y_987j). Following Han & Gould ( |1995| ) the disper- 
sion velocities in the coordinate system given by the principal 
axes of the bulge ellipsoid are (ov , oy, oy) = (113.6, 77.4, 
66.3), where the mean line of sight dispersion velocity is nor- 
malized to 110 km s . The projected velocity dispersions in 
Galactic coordinates (such that x is along the line of sight, x — y 
is in the Galactic plane and z perpendicular to it) are then com- 
puted to be (a x , a y , a z ) = (110, 82.5, 66.3). In the following 
it is convenient to write the distribution function f(v y , v z ) in- 
troducing polar coordinates in the v y — v z plane. This way in 
polar coordinates the transverse velocity distribution turns out 
to be 



fa B (y,<&)dvdd = - e~ u dvdti . 

TT 

The transverse velocity is then given by 



v T = vJ2{alcos 2 d + ofshri?) . 



(11) 



(12) 



d is the polar angle and the variable v is dimensionless. 

In the halo we consider a Maxwellian velocity distribution 
with typical dispersion velocity of vh = V2cth — 220 km s _1 



(Paczyhski 1991). The transverse velocity distribution is then 
given by 



fa,H(v T )dv T 
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For lenses belonging to the disk, the velocities of the ob- 
server and the source transverse to the line of sight are vo 
and v s , respectively, and the transverse velocity of the mi- 
crolensing tube of radius Re (where the Einstein radius Re 



will be defined in Sect. 3) at position xD (0 
Vt(x) = (1 — x)vq + a;v s . Its absolute value is 



< x < 1) is 



vt(x) = (14) 

V(l - a;) 2 |vo| 2 + z 2 |v s | 2 + 2xll - x)|v ||v s |cos ■& , 

where v = v Q cos I and v s are the solar and the source ve- 
locities transverse to the line of sight and d the angle between 
them, whereas I denotes the Galactic longitude. The solar ve- 
locity transverse to the Sun-Galactic centre line is denoted by 
v Q , with |v©| ~ 220 kms" 1 . The distribution of the trans- 
verse velocity of the lens thus is 



gr( v T, dvxd'd 



1 



vt exp 



(15) 



" U D \ V D 

with a velocity dispersion vp ~ 30 km s _1 in the disk (Pac- 
zyhski [1991]). 



Due to the relative position 7 Nor and 1} Mus are approx- 
imately moving towards us, whereas 7 and (3 Set are moving 
away from us. Hence, we will neglect the dispersion in the 
transverse velocity for sources located in these particular fields 
towards the spiral arms. We take into accout the dispersion ve- 
locity of the lenses in the disk using an expression as given in 
Eq.(p3|) but with a velocity dispersion vp ~ 30 km s~ . 

3. Mass and time moments 

3.1. Optical depth and differential number of events 

To estimate the microlensing probability one introduces the op- 
tical depth T opt defined as 

rl 4vrG 



' opt 



-p(x)D x(l — x) dx , 



(16) 



with p(x) the mass density at the distance s = xD from the 
observer along the line of sight. r opt is the probability that a 
source is found within the Einstein radius Re of a lens, where 
Re is given by 



Ri 



AGMD 



x(l — x) = r E px(l — x) 



(17) 



with x = s/D; D and s are the distances to the source and the 
lens, respectively. 

The optical depth is independent of the mass function of 
the lensing objects. To get the number of microlensing events 
one introduces t he diff erential num ber of microlensing events 
(De Rujula et al. |l99lj Griest |l99l|) : 



dN ev 



(18) 



dn 



N^t obs UTH2v T fa(vT)DrE\/^x(l - x) ^d^dv T dx , 

dfj, 

assuming an experiment that monitors stars during a period 
t a i, s . This quantity yields the number of microlensing events 
with a magnification above a certain threshold Ath (where 

A TH = ("th + 2)/(uth04h + 4 ))> A ™ = 1-34 for 
uth = !)■ 
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3.2. Method of mass moments 

A systematic way to describe the relevant quantities in a mi- 
crolensing experiment is the method of mass moments (De 
Rujula et al. 1991), which are defined as 



dfi e„(/x) -j— p r ' 



< p" 1 > 
with m = (n + l)/2 and 



_ [dN£M e(T) 
fdN* v {p)r" 



(19) 



(20) 



The sampling efficiency e(T) is given by the experiment and 
r = (vi/rE)T, where T is the time scale of the observed event 
(for which we adopt the definition that it is the time needed 
to cross the Einstein radius) and V{ stands for the dispersion 
velocity in the bulge, disk or halo. e n (p) is called efficiency 
function and measures the fraction of the total number of mi- 
crolensing events that, for a fixed MACHO mass M — pM®, 
meet the condition T m . m ^ T ^ T max . dN* v (p) is defined as 
dN ev with, however, dno / dp assumed to be a delta function 
6(p — p). In other words, eo(/i) indicates how efficient the ex- 
periment is to detect a MACHO with a given mass M = pM & . 

< p m > is related to the cumulative n th moment of r con- 
structed from the observations as follows 



events 



(21) 



Due to the insertion of e„(/i), the theoretical expression for the 
time moment 



< t" >= 



dN ev e n (p)r r ' 



factorises as follows 

< T n >= Vu TH "f{m) < p m > 
with 

V = 2N,t obs Dr EVl . 



(22) 



(23) 



(24) 



7(m) is a quantity defined by the model, which depends on the 
spatial and the velocity distributions. 

The mass moment < p m > is related to < r" >, as given 
from the measured values of T, by 



< p >= 



< t > 
Vu TH j(m) 



(25) 



Some moments are directly related to physical quantities: the 
number of events is N°^ s =< r° >. The mean local density of 
MACHOs (number per pc 3 ) is < p° > and the average local 
mass density in MACHOs (solar masses per pc 3 ) is < p 1 >. 
Thus the mean MACHO mass is given by 



< p 1 > _ < T 1 > 7(0) 

< P° > < T- 1 > 7(1) 



(26) 



in units of solar masses, where < r 1 > and < r _1 > are deter- 
mined through the observed microlensing events. The average 
event duration < T > can be expressed as follows 



< T > = 



r E < p 1 > 7(1) 
Vi < p 1 / 2 > 7(1/2) 



(27) 



where again i stands for the different components of the Galac- 
tic model. 



3.3. Sources in the Galactic centre 

Due to the extension of the bulge one has to take into account 
that the number density of possible source stars varies with dis- 
tance D. This effect has to be taken into account in the calcula- 
tion of the optical depth and the n umber of microlensing events. 
Following Kiraga & Paczynski ( 1994 ), the volume of space 
varies with distance as D 2 dD and the number of detectable 
stars as assuming that the fraction of stars brighter then 
some luminosity L is proportional to L 13 , where (3 is a constant. 
In the following we will use a power-law luminosity function 
~ L^ 1 (i.e. (3 = —1), which is appropriate for main sequence 
source stars. An observable quantity towards the bulge, which 
we shall generic ally denote as g E , is then obtained as 



9b = 



/ °° <p(D)g(D)dD 
J™v(D)dD 



(28) 



with (p(D) = pbD 2+213 and g(D) an observable quantity as 
for instance t(D) or N ev (D). ps is the number density of stars 
as given in Eq.QlJ) and also varies as a function of D. 

Accordingly, Eq.(|22|) and the subsequent relations which 
define the moments get modified, due to the additional inte- 
gration over the distance D. It is, therefore, necessary to use 
the time T instead of r (which via r E depends on D). Other- 
wise, however, the calculation goes through similarly and the 
factorization of the mass moments is again achieved. Indeed, 
the D dependence in the efficiencies e n (p) is negligible as we 
checked numerically. 

Eq.(E2b then becomes 



< T n >= V(n)7(m) < p m > 
with 
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(30) 
(31) 



where N B = / °° tp(D)dD = / °° p B D 2 + 2 f j dD. 

On the other hand < T n > is determined from the obser- 
vations by < T n >= Y^, T n (T given in days). Thus, the mean 
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mass towards the bulge turns out to be 

< /J 1 > fee) 2 < T 1 > 7(0) 



< /i u > 



4GM < T- 1 > 7 (1) 



5.35 



Vkms- 1 / 



2 < T 1 > 7 (0) 
< T- 1 >7(1) 



and the event duration is 



<T> = 



2\/GM& 7(1) <^> 

UjC 7(1/2) < /xV2 > 

-1 7 (1) <M X > 



(32) 
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7(1/2) < M V2 > 



days . 



We now turn to the computation of the 7(771), which en- 
ter into Eq.(pl|). For microlensing events towards the bulge, 
we have to distinguish between two cases: a) both the source 
and the lens are located in the bulge (bulge-bulge event), b) 
the source belongs to the bulge and the lens to the disk (bulge- 
disk event). We suppose the velocity distribution of sources and 
lenses in the bulge to be Maxwellian as given by Eq.jl^). In this 
case the relative motion between observer and source can be 
neglected and the 7(m)-function as defined in Eq.(^) through 
the mass moments is given by 



7-BB (to) = 3(2 - m)H(m) , 
where 

H(m)= [ (x(l - x)) m ^^-dx 



P°B 
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For lenses belonging to the disk, the transverse velocity distri- 
bution is given by Eq.([l5|), and we get 

/>27T />OC 

jBD(m) = 2 / d-d dv s f T , B (v > s ,fi) 
Jo Jo 

Pd{x) 



[x(l - x)) m e"" 



/ dx 

10 P°D 
poo 

x / dyy z - 2m h{2 W )e-y 2 
Jo 



(37) 



with y = vt/vd, and 7/ = v t /vD- The index D labels the disk, 
vd ~ 30 km s _1 is the dispersion velocity in the disk and Iq is 
the modified Bessel function of order defined as 



Io(x) 



1 



d(3 



(38) 



The source stars belong to the bulge and are supposed to follow 
the distribution /t,b(i> s , $) as given in Eq.dTlh. 



3.4. Sources in the spiral arms and in the halo 

If the source as well as the lens is located in the disk, the small 
relative transverse velocity between them can be neglected, as 



explained above and we can assume a Maxwellian velocity dis- 
tribution as given in Eq. (|l3]). We can perform the ■& integra- 
tion and the corresponding 7 (to) -function turns out to be as in 
Eq.rfPb, with Eq.fsel) replaced by 



H(2 



m)= / 

Jo \ V D 



fT,H{vr)dvT 



T(2 



TO 



(39) 



where vu — 30 km s _1 is the dispersion velocity in the disk. 
Sources in the halo, like LMC, SMC or dwarf galaxies are 
more distant, and thus the relative motion between observer 
and source can also be neglected in leading order. Assuming in 
addition a Maxwellian velocity distribution the resulting 7(771.)- 
function turns out to be as in Eqs. (|^) and (^), with re- 
placed by Vh, the dispersion velocity in the halo. 

4. Lensing towards the Galactic bulge and some spiral arm 
regions 

Let us apply the formalism presented above to different lines 
of sight towards the bulge and some fields in the spiral arms. 
In the bulge we focus on the direction of Baade's Window (I = 
1°, b = —3.9°) and in the spiral arms towards 7 Nor (I = 
—28.8°, b = —2.7°), where some microlensing events have 
already been observed. 

4.1. Optical depth 

In Fig . [TJ the optical depth r is shown as a function of the Galac- 
tic longitude I for a fixed value of the latitude b = —3°. The 
approximate position of 7 Nor, 7 and (3 Set and Baade's Win- 
dow are indicated. The idea is just to give an impression of the 
global behaviour of the optical depth towards Baade's Window 
and the spiral arms. The values towards these targets are given 
below. 

The dashed line tb b is the contribution due to lenses and 
sources located in the bulge for which we have taken into ac- 
count the varying distances. One clearly sees that it is the major 
contribution to the optical depth in the range —8° < I < 10°. 
The tilt angle $0 (for which we assume <E>o = 20°) of the bar 
major axis with respect to the Sun-Galactic centre line leads to- 
wards the bulge to a slight asymmetry in the optical depth with 
respect to I = 0°. 

The curve labeled by tb is the total optical depth, which 
contains also the contribution due to lenses located in the disk 
and the source star in the bulge for which we again take into 
account the varying distances. 

Towards Baade's Window (I = 1°, b = -3.9°) the total 
optical depth tb = 1.9 x 10~ 6 is the sum of the disk tb = 
0.7 x 10~ 6 and the bulge t B b = 1.2 x 10~ 6 contributions. 
The bulge contribution is about 1 .7 times more important than 
the disk. This value is in reasonable agreement with the optical 
depth found b y the M ACHO group: t b = (2.4 ± 0.5) x 10~ 6 
(Alcock et al. 1997a). The OGLE collaboration found a some- 
what higher value: t b = (3.3 ± 1.2) x 10~ 6 (Udalsky et al. 
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-30.0 -20.0 -10.0 0.0 10.0 20.0 30.0 

Longitude 1 

Fig. 1. The optical depth r (in units of 10 -6 ) for the different com- 
ponents of the Galaxy as a function of the Galactic longitude I for a 
fixed latitude b — —3°. For the disk contribution (td) the distance 
is fixed at 8 kpc, while for the bulge contribution (r_e) we integrated 
over the distance distribution of the lenses and sources, tbb denotes 
the fraction of r for source and lens located in the bulge (where the tilt 
angle $o is taken to be 20°). The dotted lines are just an interpolation 
between the bulge and the disk regions. The approximate positions of 
7 Nor, 7 and ft Set and Baade's Window (BW) are indicated. 



1994). We notice that it has been pointed out that the contri- 



bution of unresolved stars might be quite significant and as a 
result i mply that the measured optical depth is overestimated 
(Alard 1997 ). If this is the case the measured value has to be 
regarded as an upper limit to the true value. 

By varying the tilt angle from 15° to 30° the optical depth 
T BB for the bulge decreases from 1.33 x 1(T 6 to 0.97 x 1(T 6 
(and the total optical depth from 2.02 x 10~ 6 to 1.68 x 10~ 6 ). 

Outside the bulge region both the source and the lens be- 
long to the disk population and the total optical depth is la- 
belled in Fig. [l]by 775. In this case we assume all source stars 
to be located at the same distance of 8 kpc. 

Observations towards 27 different fields in the spiral arms 
are included in the scientific program of the EROS II (Derue 



et al. 1999) and the MACHO teams. The 27 fields observed by 



EROS II belong to 4 dense regions of the spiral arms: 7 Set, /3 
Set, 7 Nor and -d Mus. The star density lies between 400000 
and 700 000 stars per field, leading to about 10 millions observ- 



able stars (Mansoux 1997 ). Their average Galactic coordinates, 



distances, observed number of stars, number of events found so 
far during an observation time of 1 .7 years and the calculated 
optical depth are reported in Table [|. 

EROS II found altogether 3 events (see Table [j]) and esti- 
mated an optical depth, averaged over the four directions, of 
0.38i°;ig x 10~ 6 (Derue et al. |1999| ). For comparison we 



4.2. The mass functions 

Towards the bulge we compute the mass functions using the 
first year data of the MACHO collaboration taken in 1994. Ob- 
serving during 190 days star fields with a total of 12.6 mil- 
lion source stars towards Baade's Window, they found 40 mi- 
crolensing events (neglecting the double lens event). 

Due to the limited amount of data, we make the Ansatz 
dno/dfi ~ C}{ai)n~ ai and use then the mass moments to 
determine the parameters a and /i m , n . C}{ai) is defined by 
Eq.(^|) and we assume a maximal mass (fx max ) for the stars of 
~ 10M@. The upper limit fj, up of the integration in Eq.(|l9|) is 
set equal to 1 M®, since more massive and thus brighter stars 
will not contribute much to microlensing events. To evaluate 
Eq.(|l9|) we first have to calculate the efficiency functions de- 
fined by Eq.(pO|) using the sampling efficiency e(T) as given by 
the MACHO team for their first year events towards the Galac- 
tic bulge. 

We notice that one should also take into account the con- 
tribution due to unresolved stars, which for the bulge might be 
relevant and induces changes in the optical depth and the mi- 
crolensing rate (Alard 1997; Han 1997] ). Indeed, if a faint unre- 
solved star is lensed close enough to a resolved star, the event 
will be seen by the microlensing experiment and attributed 
to the brighter star. The blending biases the event towards a 
shorter duration, leading to an overestimate of the amount of 
these events. As a result the mass function will be shifted to- 
wards lower values of /i m ;„. A quantitative estimate of the in- 
duced error on our results is, however, difficult, since it would 
require a detailed knowledge of the data analysis and the deter- 
mination of the sampling efficiency. Hopefully, this informa- 
tion will be available when the new microlensing observations 
towards the bulge will be made available by the observational 
teams. 

The mass moments as computed by Eq.(|l9|) are compared 
to the corresponding time moments derived from the observed 
microlensing events through the expression given in Eq.(|25|). 
Thus with the mass moments, < fjP > and < /i 1 / 2 >, we get 
two equations, which then are solved with respect to the un- 
known quantities a and fj, m in- Further moments can be used for 
instance to check whether the values we obtain for a and £i m j„ 
are consistent. Since we do not know the fraction of the ob- 
served events due to lenses located in the bulge, respectively in 
the disk, we assume that this fraction is equal to the correspond- 
ing ratio of the optical depths hence tbb/tb — Ns/Ntoti— 
25.1/40). This assumption is well fulfilled if the mean event 
duration of the bulge and disk events are the same. We will then 
a posteriori verify if this condition is matched by the results we 
get and thus see whether our assumption is consistent. 

With the above procedure we get for the parameters of the 
power-law mass function of the lenses in the bulge the values 
a = 2.0 and fj, min = 0.012 with C B = 0.34 Af Q pc" 3 . Simi- 
larly we can use the disk events to determine the mass function 
of the disk population to find again a = 2.0 and instead a min- 



get an averaged optical depth of r = 0.42 x 10 using our imal mass of 0.021 M Q with C D = 9.86 x 10 M pc 



model, which is in good agreement with the EROS II value. The average durations we find are (Tb) 



20.5 days and 



7 



(To) = 19 days, respectively. Therefore we see that the as- 
sumption (Tg) ~ (Td) is quite well verified. 

Varying the tilt angle $ between 15° and 30° the slope a 
for the bulge mass function slightly changes from 1.9 to 2.15 
and the minimal mass fi m i n from 0.009 to 0.023 whereas the 
number of bulge events decreases from 25 down to 22.7. The 
mean event duration for the bulge (Tb) changes accordingly 
from 22. 1 to 19.1 days. Similar small changes occur in the cor- 
responding numbers of the disk events. We thus see that our 
results are robust for reasonable values of the tilt angle. 

Due to the above considerations on the unresolved stars, 
the values for /Lt m j„ we find have to be considered as a lower 
limit on the true values. We notice that a change in n m i n af- 
fects also the value of a. Similar studies have also been per- 
formed with different methods by Han & Gould (1995) and 
Peale (1998,1999). They find somewhat different mass func- 
tions for the bulge, but again given the small amount of data at 
disposal the errors are still large and so no definite conclusions 



can be drawn. Recently, Reid et al. (1999) using observations 



from the Deep Near-Infrared Survey (DENIS) and the 2 Mi- 
cron All-Sky Survey (2MASS) have suggested the presence of 
a local unknown population of free brown dwarfs with a mass 
function with minimal value as low as 0.01 Mq. If confirmed it 
may well be plausible to find this population also in the bulge. 

4.3. Mass and time moments towards the bulge 

For the calculation of the mass and time moments towards the 
bulge we have to distinguish between the two cases, whether 
the lens is located in the bulge or in the disk, and to choose 
accordingly th e ap propriate velocity distribution as given by 
Eq.(|TT|) or Eq.(|r5|). The number of events and the event dura- 
tion are defined by N°^ s —< T° > and Eq.(^4| ), respectively, 
where the moments of T are given by Eq.(29). 

Figure ^ shows the number of expected microlensing events 
as a function of the longitude for a delta mass distribution with 
j«o = 1, assuming to monitor 10 6 stars during an observation 
time of 1 year. The efficiency is taken to be 1 . Again, as in the 
plot of the optical depth, it shows the global behaviour of the 
expected number of events towards the bulge and the spiral arm 
regions. The distance towards the spiral arms is kept fixed at 8 
kpc and the latitude at b = — 3°. 

N eVt B is the total number including also the contribution 
due to lenses located in the disk and sources in the bulge. For 
the latter we take into account the dependence on the distance 
as illustrated in Eq.(27). For the number of events towards the 
spiral arms, N eV: n, the source is located in the spiral arms and 
the lenses belong to the disk population. 

In Fig. 4 we plot the efficiency functions which are derived 
from the sampling efficiency e(T) as given by the MACHO 
team towards Baade's Window. We recall that for instance N ev 
scales like ~ v fi" 1 ^ 2 , whereas < T > as v^ 1 (i 1 / 2 . N eVt BB 
is the contribution due to the bulge (i.e. events with source and 
lens in the bulge). 

Figure || shows the expected number of microlensing events 
assuming the power-law mass functions we found with the pa- 
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Fig. 2. The number of events for a delta mass distribution with /io = 1 
as a function of the longitude with fixed latitude b — —3°. The total 
exposure is assumed to be 10 6 star-years. The N ev> BB curve (source 
and lens in the bulge) reflects the position of the bulge (tilted by $o = 
20°), which is important in the range —8° < I < 10°. In this range, 
the total number (denoted by N cv ,b) is dominated by the events with 
the lens located in the bulge. Outside of this range the events N eVt o 
are due to sources located in the spiral arms and lenses in the disk 
population. The efficiency function eo(m) ls assumed to be 1 in this 
plot. 



rameters as mentioned in Sect. 4.2. For this plot the same ef- 
ficiency e(T) as given by the MACHO team towards Baade's 
Window is assumed for all directions. All the other parameters 
are the same as the ones used in Fig. |[ 

Using the values for the various parameters as given in Sect. 
2 and power-law mass functions with a = 2.0, fi m in,B = 
0.012 for the bulge and [i m in.D — 0.021 for the disk popula- 
tion, respectively, the number of events towards Baade's Win- 
dow are consistently reproduced to be N ev = 40.5. This is just 
the sum of N ev> BB — 25.0 and N eV: r> = 15.5, which we 
get when considering 12.6 million source stars observed dur- 
ing 190 days as in the MACHO first year data. The mean event 
duration, as given by Eq.(j34]), is < T >b= 20.5 days for bulge 
events and < T >d~ 19.0 days for disk events. These dura- 
tions have to be compared with the actual observations of the 
reported 40 events with < T f, s >= 19.9 days. We thus see 
that the values for the durations we get using our model are in 
good agreement with the observed one. For the mean mass we 
get 0.091 M Q for the lenses in the bulge and 0.124 M Q for the 
lenses in the disk, respectively. Due to the sampling efficiency 
this value is actually an upper limit of the true value (when 
not taking into account the problem of the unresolved sources). 
Just for comparison the theoretical value we would get for the 
average mass in the bulge assuming the power-law mass func- 
tion as mentioned in Sect. 4.2, with e = 1, is 0.054 Mq. As 
expected, this value is somewhat lower than the one inferred 
from the data. 
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< I > 


< b > 


D 
kpc 


AT, 
(xlO 6 ) 


N oe 


T 

(days) 


Ttot 

(xlO 6 ) 


N ev 


< T > 
(days) 


7(0) 


7(1/2) 


7(1) 


7 Set 


18.6 


-2.6 


6.5 


1.70 


1 


73 


0.41 


0.60 


76.9 


1.62 


0.563 


0.269 


/3 Set 


27 


-2.5 


6.5 


1.96 







0.38 


0.63 


77.4 


1.47 


0.516 


0.248 


7 Nor 


-28.8 


-2.7 


8.0 


3.01 


2 


98&70 


0.56 


1.29 


86.5 


1.42 


0.504 


0.244 


# Mus 


-53.6 


-1.8 


7.0? 


1.77 







0.32 


0.47 


80.9 


1.07 


0.379 


0.184 



Table 1. The distance D, the num ber of observed stars TV, and the number of observed events N oe with their duration T are given for 7 Set, f3 
Set, 7 Nor and 1? Mus (Derue et al. 1999). The optical depth Ttot, the number of events N ev and the mean event duration < T > are computed 
adopting the model outlined in Sect. 2 and a power-law mass distribution with a = 2.0 and Umin = 0.021. The disk dispersion velocity is 
vd = 30 km s _1 . We assume to observe N, stars during an observation time of 1.7 years with a sampling efficiency as given by the EROS II 
team. 





Ttot 

(xlO 6 ) 


N ev 


< T > 
(days) 


A 4 


7(0) 


7(1/2) 


7(1) 


Bulge 


1.20 


25.0 


20.5 


0.091 


3.82 x 10 2 


1.59 x 10 4 


9.58 x 10 5 


Disk 


0.70 


15.5 


19.0 


0.124 


2.83 x 10 5 


2.02 x 10 6 


1.94 x 10 7 


Total 


1.90 


40.5 


19.9 











Table 2. The optical depth, the number of events and the mean event duration towards Baade's Window are computed following the model 
outlined in the text. For the computation of N ev and < T > we inserted the computed power-law mass functions with a = 2.0 and fj, m in 
(),(! I ? for the bulge and ji m in = 0.021 for the disk. We assume an observation time of 0.52 years and 12.6 x 10 6 stars. 
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Fig. 3. The number of events for power-law mass functions with 
a = 2.0 and i± m i n ,B = 0.012 for the bulge and fx m in,,D = 0.021 
for the disk population, respectively. The efficiency function of the 
MACHO collaboration is taken into account in the calculation to be 
able to directly compare with the experiment. The total exposure is 
taken to be 10 6 star-years. 



4.4. Mass and time moments towards the spiral arms 

The corresponding results for the various directions in the 
spiral arms are summarized in Table [j] To calculate the 
number of expected events we assume the power-law mass 
function obtained for the disk and the number of observed stars 
(indicated in Table [jj) and the observation time of 1.7 years as 



0.50 




M 

Fig. 4. Efficiency functions e n (/x) (n = —1,0, 1) derived from the 
sampling efficiency towards Baade's Window of the MACHO collab- 
oration. 

given by the EROS II experiment. Moreover, we adopted the 
sampling efficiencies quoted by EROS II (Derue et al. 1999). 

For the 7 Nor direction we studied how the results change 
by varying the various parameters, which enter the calculation 
of the number of events and their duration. As a first point we 
notice that changing the maximal mass fi ma x of the source stars 
or the maximal mass /j, up of the lenses in the power-law model 
has a very minor influence on the subsequent results and can, 
therefore, be neglected. By varying the minimal mass (i m in be- 
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tween 0.005 and 0.05 keeping all other parameters fixed, N ev 
decreases from 1.7 to 1.0 and the mean event duration < T > 
increases from 63 to 106 days. If the slope value a of the 
power-law varies between 1.7 and 2.5 (for fJ, m i n = 0.021), 
the number of events increases from 0.95 to 1 .9 events with a 
corresponding duration decrease from 101 to 69 days. 

If we assume instead of a power-law a delta mass distri- 
bution with most probable masses //on = 0.124 (or /j,qb = 
0.091) we get 3.1 (3.5) events with a duration of 87 (76) days. 

Another important parameter is the dispersion velocity in 
the disk. Changing its value between 20 and 40 km s _1 , N ev 
increases from 0.86 to 1.7 and < T > decreases from 130 to 
65 days, respectively. 

We thus see that the small number of observed events does 
not yet constrain very much the different parameters entering 
the model. Nevertheless, once more data will be available the 
allowed range for the parameters will be narrowed considerably 
and in this respect the method of mass moments will be very 
useful leading to a much better knowledge of the structure of 
the bulge and the disk of our Galaxy. 



0.40 



0.30 



0.20 



o.io 



0.00 




Fig. 5. Efficiency functions e n (/x) (n = —1,0, 1) derived from the 
sampling efficiency towards the LMC of the MACHO collaboration. 



5. Lensing towards halo targets 

Two important unknowns in Galactic models are the shape of 
the halo and the velocity distribution of the halo objects. Al- 
though it is widely believed that dark matter halos are not 
spherical, a systematic explanation how to generate deviations 
from spherical symmetry is not yet available. To investigate 
how the presence of massive Galactic companions such as the 
Magellanic Clouds may induce a flattening of the halo, and to 
test how such a presence could influe nce th e isotropy of the 
dispersion velocity (Holder & Widrow, 1996 ) we performed a 
N-body simulation which yields, as a result, the density and ve- 
locity distribution of the halo objects. For the details we refer 
to the Appendix. Of course our results have to be taken as an il- 
lustration, given also the crude approximations we use. For the 
targets towards different halo directions we will adopt the flat- 
tening parameter qn and the dispersion velocities as obtained 
from the numerical simulation. These values are tabulated in 
Table U 

5.1. Magellanic Clouds 

The optical depth towards the LMC for the spherical halo 
model is t lmc - 5.33 x 10~ 7 (4.93 x 10~ 7 due to the 
halo and 0.40 x 10 -7 due to the disk contribution), where 
we used the standard parameters for the LMC (I = 280.5°, 
b = —32.9°, D — 50 kpc). For the halo core-radius we adopt 
Rc = 5.6 kpc. The measured value, as reported by the MA- 
CHO collaboration, is: tmacho — 

x 10~ 7 , which cor- 
responds to about 50% of the above predicted value for a stan- 
dard spherical halo (Alcock et al. 1997ai. Similarly, towards 



The MACHO team found two microlensing events towards 
the SMC, one being a binary event (Alcock et al. |l997c] Pa- 



the SMC (2=302.8, 6=-44.3, D=63kpc) the calculated optical 
depth is tsmc - 7.39 x 10~ 7 (7.08 x 10~ 7 and 0.31 x 10~ 7 
for the halo and disk contribution, respectively). 



lanque-Delabrouille et al. 1997, Udalski et al. 1997). Using a 
simple estimate for the total exposure and assuming for e the 
corresponding LMC efficiency, they estimate the optical depth 
towards the SMC to be roughly equal to the optical depth to- 
wards the LMC. A similar conclusion has been also reached 
by the EROS II team, which quotes a value r ~ 3.3 x 10 -7 
(Perdereau 1998). To get a ratio tlmc/ t smc of about 1, a 
halo flattening of almost qn = 0.5 is required. However, it 
might well be that the events found so far are due to lenses 
in the SMC itself (in particular this is the case for the binary 
event), in which case the above conclusions on r are no longer 
valid. 

From the 6 events (excluding the binary lens event and 
an event which is considered being only marginally consistent 
with microlensing) published by the MACHO team for their 
first two years observations, we find with the method of mass 
moments f or a sp herical halo model an average mass of 0.26 
Mq (Jetzer 1996| ). If we assume for the dispersion velocity 191 
km s _1 as found in the N-body simulation instead of 210 km 
s _1 as used for the spherical halo, we get for the average mass 
a slightly lower value of 0.22 Mq. 

Adopting the values for the first two years observation of 
the MACHO team towards the LMC, namely a total exposure 
of 2.1 years for 8.5 x 10 6 stars with the detection threshold 
Uth = 0.661 and the calculated efficiency functions e„(//) 
from the MACHO sampling efficiency, we compute the ex- 
pected number of events and the mean event duration for a 
qn = 0.8 flattened halo. A delta mass function with fijj = 0.22 
is used for the halo and the computed power-law mass function 
{p, min = 0.021, a = 2.0) for the disk. We find N ev = 12.2 
(where the halo contributes with Njj = 10.4 and the disk with 
N D = 1.8) and < T >lmc= 32.5 days « T > H = 34.1 
days and < T >d= 23.1 days, respectively). Using a spher- 



10 



Microlensing Events and Dark Objects by Mass and Time Moments 



ical halo with vh = 210 km s _1 we find N ev — 14.5 and 
< T >lmc= 30.3 days. We see that the flattening of the halo 
and the change in the velocity dispersion reduces, although 
slightly, the number of expected events, while the average du- 
ration increases. 

The value for N ev is valid assuming a halo made entirely 
of MACHOs. For a 50% contribution of MACHOs as implied 
from the optical depth we thus expect about 6 events, which 
compares well with the 6 observed events. The event duration 
has to be compared with the measured 35.5 days of the 6 ob- 
served events. It is noticeable that we expect about two events 
as due to the disk contribution. In addition LM C self- lensing 
may also contribute with some events (Salati et al. 1999 ), which 
we do not take into account in our calculation. Indeed, the frac- 
tion of the optical depth due to self-lensing is st ill co ntrover- 
sial, although a recent analysis of Gyuk et al. ( 1999 ) comes 
to the conclusion that it contributes at most 10-20% to the ob- 
served optical depth. Thus the bulk of the events should still 
come from MACHOs located in the halo. Given this mixing of 
different populations in the events one has to be careful with 
the interpretation of the average mass, since it is derived under 
the assumption that all events are due to MACHOs in the halo, 
which is probably not the case. 

5.2. Dwarf galaxies 

The LMC and the SMC are the only two targets in the halo with 
a sufficiently high number of stars which can be resolved such 
that the method of microlensing can be applied. Other targets 
such as the satellite dwarf galaxies, can still be used for mi- 
crolensing observations however, due to their size and distance 
much less stars can be resolved. One could envisage in this case 
to use the pixellensing method instead. With this method one 
can look at rich fields of stars which are located further out 
and thus are not resolvable. The method, proposed by Crotts 
\\992\ and Baillon et al. ( |1993| ), has proven to work with M31 
and so it is conceavable to extend the observations to some of 
the several satellite dwarf galaxies of the Milky Way. Such ob- 
servations would be very useful, since they would allow to test 
other directions as the ones towards the LMC and the SMC. 

In Table || we report the relevant quantities for some pos- 
sible targets in the halo and the neighbourhood of the Milky 
Way, for whi ch we assumed an extension of the halo of 150 
kpc (Bahcall 1996 ) and a flattening parameter of qn = 0.8. 
For the number of events and the mean event duration a delta 
distribution with the mass value //o = 1 is assumed. The disper- 
sion velocity is taken from the simulation presented in the Ap- 
pendix. The total exposure is taken to be 10 6 star-years and the 
efficiency is set equal to 1 . The number of events is determined 
under the assumption of a halo made entirely of MACHOs. 

In comparison for a standard spherical halo with vh = 210 
km s _1 we would expect about 20% more events, but with a du- 
ration of about 7% shorter. Similar results have also been found 



by De Paolis et al. (1996) in the framework of halo models with 



anisotropy in the velocity space. 



Other possible targets from which interesting Galactic 
structure can be extracted are globular clusters, for which we 
refer to the papers by Gyuk & Holder (1998), Rhoads & Mal- 
hotra (1998) and Jetzer et al. (1998). 



6. Discussion 

We studied in detail microlensing towards different lines of 
sight which are promising to obtain information about the 
structure of our Galaxy using the method of mass moments. 
With the first year MACHO data in the direction of the Galactic 
centre we determined the mass functions of lenses in the bulge 
and the disk assuming a power-law. The obtained mass func- 
tions are both slightly less steep continuations of the Salpeter 
IMF down to the brown dwarf region, with a minimal mass of 
~ 0.012 M Q for the bulge and ~ 0.021 M for the disk with 
a slope a ~ 2.0 for both cases. By varying the tilt angle of the 
bulge in the range from 15° to 30°, which is suggested by the 
most recent models, the inferred mass functions do practically 
not change. 

As next we computed the expected number of events and 
their duration for different targets in the spiral arms, which 
have been explored by the EROS II collaboration. Assuming 
the same mass function as derived for the disk we get event 
durations which are in good agreement with the values found 
for the 3 events observed towards 7 Nor and 7 Set. On the other 
hand, the calculated number of events turns out to be somewhat 
smaller, although given the few events found so far the uncer- 
tainties are still large and, moreover, the theoretical values vary 
in a significant way even for slight changes of the parameters 
such as the velocity dispersion. 

Microlensing towards fields in the spiral arms is an impor- 
tant tool to explore the structure of the disk and especially its 
mass function and dispersion velocity. The method of mass mo- 
ments can be used to get in a systematic way important infor- 
mation from the data, as soon as a sufficient number of events 
will be observed, which is certainly the case towards Baade's 
Window. It is also important to further develop the Galactic 
models, especially the modelling of the spiral arms and the mo- 
tion of the stars inside them. Comparison then with the data will 
lead to stringent constraints on the model. 

To obtain a qualitative understanding how the presence of 
massive companions may influence the Galactic halo, we per- 
formed a N-body simulation of the halo from which we ob- 
tained a halo flattening parameter qn — 0.8 as well as an idea 
how the dispersion velocity can vary as a function of the ob- 
servational direction. The parameters from the simulation were 
then used to compute the optical depth, the expected number 
of events and the duration for several targets in the halo of our 
Galaxy. 

Although most of the results we obtained have to be con- 
sidered as preliminary and will soon improve due to many new 
events which will be available, we find that within the adopted 
model for the Galaxy the agreement between the computed val- 
ues for the optical depth, the expected number of microlensing 
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target 


distance 


1 


b 


Ttot 


disp. vel 


N ev 


< T > 


7(0) 


7(1/2) 


7(1) 




(kpc) 






(xlO 7 ) 


(kmr 1 ) 




(days) 


(xlO 2 ) 


(xlO 2 ) 


(xlO 2 ) 


LMC 


50 


280.5 


-32.9 


4.8 


191 


1.5 


73 


28.48 


8.60 


3.71 


SMC 


63 


302.8 


-44.3 


6.1 


195 


1.8 


78 


25.69 


7.34 


3.04 


NGC 205 


690 


120.7 


-21.1 


7.3 


344 


2.8 


59 


6.82 


1.64 


0.62 


NGC 147 


690 


119.8 


-14.3 


7.8 


195 


1.9 


96 


8.02 


1.81 


0.65 


Fornax 


230 


237.3 


-65.7 


6.2 


195 


1.4 


102 


7.72 


1.65 


0.57 


Draco 


60 


86.4 


34.7 


5.1 


192 


1.5 


77 


22.73 


6.54 


2.74 


Ursa Minor 


80 


104.9 


44.9 


5.0 


184 


1.3 


89 


15.06 


3.95 


1.57 


Ursa Major 


120 


202.3 


71.8 


5.3 


175 


1.2 


106 


9.21 


2.10 


0.76 


Leo I 


230 


226.0 


49.1 


6.2 


174 


1.3 


115 


7.38 


1.60 


0.56 


Leo II 


230 


220.1 


67.2 


6.0 


176 


1.2 


113 


7.48 


1.59 


0.55 


NGC 2419 


60 


180.4 


25.3 


3.4 


202 


1.1 


70 


15.53 


4.33 


1.80 


Sagittarius 


24 


5 


-15 


12.6 


205 


5.2 


56 


143.29 


53.74 


26.92 



Table 3. Optical depth, number of events and event duration for different targets in the halo and the neighbourhood of the Milky Way. For the 
number of events and the mean event duration a delta distribution with the mass value no = 1 is assumed. The total exposure is taken to be 10 6 
star-years. The halo is assumed to be flattened with qn = 0.8 and the values for the dispersion velocity are taken from the N-body simulation. 



events and the observation is quite good, given also the uncer- 
tainties in various important parameters. 
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Appendix A: N-body Simulation of the Galactic Halo 

In this Appendix we present a N-body simulation of the galac- 
tic halo. The simulation aims to give a qualitative understand- 
ing of how and to what extent the presence of massive bodies 
like the LMC and SMC influence the shape and the velocity 
dispersion in the halo. These latter quantities are relevant for 
microlensing. We stress that our simulation has to be consid- 
ered as a first step and that the results we find should be re- 
garded as illustrative. 

A.l. The Model 

Our model consists of three massive point particles, which rep- 
resent the Galaxy, the LMC and the SMC, respectively. The 
central mass is assumed to be 5 x 10 11 M©, while we set for the 
LMC-mass 6 x 10 9 M 0) and for the SMC-mass 1.5 x 10 9 M m 



(Westerlund 1990). In the field of these three gravitating bodies 
we let evolve a halo consisting of 2.8 x 10 7 test-particles for 
about 9.5 x 10 9 years. The time resolution of the simulation is 
30 000 time steps of 10 13 s (3.17 x 10 5 years) each. The halo 
is discretised into 50 3 equidistant bins with respect to spherical 
coordinates (r, ip), hence, yielding an average population 
of ~ 200 MACHOs per bin. To keep the problem clean and 
simple, we do not yet include an extended (oblate or prolate) 
massive halo into the gravitational potential hence, we ignore 
any possible backreaction of the halo on the clouds. As a conse- 



quence of our test-particle approach all the results will be valid 
only up to a normalization constant depending on the total mass 
of the halo. 

The initial conditions for the Magellanic Clouds are chosen 
such that at the end of the simulation the locations and the ve- 
locities are close to the observed values of the LMC and SMC 



(Lin, Jones & Klemola 1995). However, since the many-body 
problem depends crucially on the initial values and we made 
some compromises to avoid pathological orbits due to close 
encounters between the LMC and the SMC, the final location 
does not exactly match the real one. 

The initial spatial distribution of the halo objects mimics an 
isothermal sphere with a cut-off radius of 100 kpc. The initial 
velocities were chosen to be Keplerian in magnitude and ran- 
dom in direction. We assume a total zero angular momentum 
halo in order to not overestimate the halo flattening due to the 
presence of the Magellanic Clouds. 

To test the reliability of the results we performed some con- 
trol runs, varying the initial conditions as well as the duration 
and the time resolution of the simulation. We found that the re- 
sults are stable unless significant anisotropics are introduced in 
the initial data. A time duration slightly less than a Hubble time 
is sufficient for the developement of a steady configuration. 



A.2. The Results 

The simulation yields the entire phase-space of all halo ob- 
jects. Hence, we are able to extract all halo parameters enter- 
ing in gravitational microlensing, which are independent of the 
sources in relative units with the exception of the lens mass 
function. 
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Fig. A.l. Isocontour lines of the MACHO-density in cylindrical 
coordinates. The final configuration is sligthly flattened with 
qn ~ 0.8. Polar orbits are significantly underpopulated. 



As it can be seen in Fig. Al, the final density profile is 
rather well described by an isothermal spheroid 



p(R, z) = 



PoH 



R 2 + z 2 /q 2 H 



(A.l) 



with qn — 0.8. An anisotropy in the disk plane is not seen, 
which is not surprising due to the model assumptions. The den- 
sity profile is slightly dimpled at the poles which again is a con- 
sequence of the LMC and SMC orbits. The radially integrated 
MACHO-density decreases towards the poles and falls below 
80% of its average value at a galactic latitude higher than ~ 70 
degrees. Although the MACHO-density also varies azimuthally 
by about 10%, a clear correlation or anti-correlation with the fi- 
nal position of the LMC or SMC is not seen. 

The angular dependence of the optical density t is quite 
similar to the radially integrated MACHO-density. Due to the 
geometrical factor when integrating along the line of sight, the 
variation is somewhat weaker. In fact, the simulation implies 
that the optical density towards the LMC should not differ more 
than ~ 20% from the average value as a consequence of the 
halo flattening induced by the presence of the LMC and SMC. 

Contrary to the density profile the velocity dispersion of the 
halo objects significantly varies with its position in the halo. 
Hence, the global distribution function in phase-space cannot 
be written as a product of a spatial and a velocity distribution. 
In fact, close to the LMC we even observe a bulk motion of 
MACHOs which form a separate halo around the LMC. To- 
wards an undisturbed direction far from the LMC or the SMC, 
the velocity distribution is rather well described by a Maxwell 
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Fig.A.2. Relative integrated MACHO-density p for an ob- 
server in the center of the halo. The distance to the source 
is set to be 50 kpc. The trajectories show the projected orbits 
of the two massive companions. At the end of the simulation 
the LMC is located at / = 80.7, b =-31.3 and the SMC at 
I = 93.7, b =-30.4. 




Fig. A.3. Relative optical density t for an observer in the cen- 
ter of the halo. The distance to the source is set to be 50 kpc. 
The trajectories show the projected orbits of the two massive 
companions. At the end of the simulation the LMC is located 
at I = 80.7, b =-31.3 and the SMC at I = 93.7, b =-30.4. 

distribution with a dispersion velocity depending on the obser- 
vational position. 

For various locations in the halo the tangential velocity dis- 
persion vh as given by the simulation is tabulated in Table |[ 
To calculate the tangential velocity dispersion, as used in grav- 
itational microlensing from the simulated spatial one, we use 

\J {v y v z cos (p sin f?) 2 + (v x v z sin ip sin #) 2 + (v x v y cos $) 2 , 

(A.2) 

where ip and i3 denote the azimuthal and the declination angle 
of the observational direction in standard spherical coordinates. 
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